EULER EQUATIONS AND TURBULENCE: 
ANALYTICAL APPROACH TO INTERMITTENCY 



A. CHESKIDOV AND R. SHVYDKOY 

Abstract. In this note we introduce in precise mathematical 
terms some of the empirical concepts used to describe intermit- 
tency in a fully developed turbulence. We will give definitions of 
the active turbulent region, volume, eddies, energy dissipation set, 
and derive rigorously some power laws of turbulence. In particular, 
the formula for the Hausdorff dimension of the energy dissipation 
set will be justified, and upper/lower bounds on the energy spec- 
trum will be obtained. 



1. Introduction 

Intermittency, if viewed as a measure of non-uniformity of the energy 
cascade, bears its signature on many statistical laws of a fully developed 
turbulence. A typical example is the power law for the energy density 
function, 

a) B( K )~4(-) 1 *, 

where e is the energy dissipation rate per unit mass, kq the integral 
scale, and d is the dimension of a set that carries turbulent energy 
dissipation. One observes similar corrections for longitudinal structure 
functions or skewness and flatness factors of velocity gradients (see the 
texts [15l[T8j[T7] for detailed discussions). Due to intermittency the de- 
terministic counterparts of the power laws appear in a form of bounds 
to account for possible corrections as in Constantin et al [21 El E] , even 
though the intermittency itself has not been mathematically defined. 
As Foias remarked [13], "The rigorous mathematical reformulation of 
the heuristic estimates provided by some widely accepted turbulence 
theories always seem to represent the upper bound of the rigorous 
mathematical estimates." At present there is no compelling analyt- 
ical evidence that the fluid equations are consistent with the full range 
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of intermittent dimensions < d < 3 (or with d < 0). On the other 
hand, experimental evidence points only to moderate deviations from 
the classical Kolmogorov theory. Based on several subjective factors, 
such as the overwhelming success of widely used hot wire anemometry, 
or high incidence rates of linear cross-sections with (2+)D fractals (see 
[16]). one is lead to believe that the values of d above 2 are more likely to 
occur in real turbulence. For Leray-Hopf solutions of the Navier-Stokes 
equation this question is directly related to the regularity problem for 
the following reason. If d is defined as a saturation level of Bernstein's 
inequality at the upper end of the inertial range (in a proper sense) 
then, as shown in [3J, a solution with dimension d above 3/2 is auto- 
matically regular. We note that such an interpretation of intermittency 
also appeared in a derivation of the dyadic [I] and continuous [5] models 
of the energy cascade. 

Without implementation of viscous dissipation the statistical theory 
offers examples of models where the dimension is introduced 
efaction exponent of Richardson's cascade (see [TJ1 Q21 HE])- Let us 
briefly recall the 0- model of Frisch, Sulem, and Nelkin [2]. Assume 
that a fluid is in a state of turbulent motion confined to the periodic 
box T| = [0,L] 3 and observed on a long interval of time [0,T]. On 
each stage of the cascade the kinetic energy is carried through by a 
set of active eddies {t q }. Let the eddied fill a volume V q . If one as- 
sumes that V q is a fraction of the preceding volume, i.e. V q ~ pV g -i, 
for some < j5 < 1, then going up the scales, V q ~ (3 q L 3 . Letting 
f) = — log 2 /3, we have V q ~ 2~ qi) L 3 and f) > 0. A heuristic inci- 
dence argument shows that the ultimate accumulation of the cascade 
occurs on a set that crosses infinitely many generations of eddies, i.e. 
A = limsupg^^ UCg, and the dimension of the accumulant A does not 
exceed d = 3 — t). A simple dimensional analysis recovers the energy 
spectrum law ([1]). Further hypothesis of a statistical uni-scaling of ve- 
locity displacements 8u(£) on A made the model unsuitable to capture 
experimentally observed exponents for higher order structure functions. 
This lead to subsequent multi-fractal ramifications of Frisch and Parisi 
(see [2HIM5]). 

Our objective in this paper is to define the above mentioned empirical 
concepts for deterministic fields, and use them to justify intermittency 
corrections to statistical laws without the use of viscous dissipation. We 
outline our basic definitions and results. Let u G L 2 ([0,T] x T|) be a 
divergence free time dependent vector field. We will use the Littlewood- 
Paley decomposition of u given by u = J2 g >~i u q ( see Section I2TTT) . We 
interpret u q as the collective velocity field of all eddies z q of dyadic size 
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£ q = L/2 q . By bridging relations between the Euleran and Lagrangian 
description of the energy flux through scales £ q we define the active 
volume as 

(9\ V - 7-3 (Kl 2 ) 3 

Here and throughout the bracket (■) denotes the average over the space- 
time domain Qt = [0, T\ x T|. The co-dimension f) is then introduced 
as the following exponential type of the sequence {V^}: 

(3) „ = , iminf igg^m 

q^oo q 

Note that f) or the dimension d — 3 — I) is a number "attached" to 
the field u without any assumptions on u. By the Holder inequality we 
have one obvious bound V q < L 3 , and in the stationary case Bernstein's 
inequality implies V q > c£ q , where c is a-dimensional. The latter implies 
d > 0, which is a natural bound. However, in the time-dependent case, 
V q can scale down faster due to possible temporal rarefaction of active 
regions giving negative values of the dimension d (see also a discussion 
in Frisch [15]). 

To define active regions that fill volume V q we find that for an eddie 
to qualify as "active" it has to have speed at least proportional to 

(M 3 > 



s n 



9 (Kl 2 )' 

The speeds will be interpreted as magnitudes of coefficients in the 
atomic decomposition of u q , and the supports of atoms will be viewed 
as eddies. The collection of those eddies will constitute what we call 
an "active" region A q . The number of active eddies is represented by 
V q /£ q . Our Theorem 13.21 shows that the Hausdorff dimension of 

A = lim sup A q 

q— >oo 

is not to exceed d as predicted by the model, and A q supports most 
of the "active" part of the field u q in terms of L 3 -averages. We will 
recover the energy law ([T]) in the form of an upper and lower bounds. 
The scaling of the second order structure function is as predicted by 
the /3-model 



(4) S 2 (£) < 53 



2 „2 
3 



d c 4 
3~ d / P\3 



+ c e[ - L 



for all £ < L, and any 5 > 0. Here and throughout, < denotes inequality 
up to a universal adimensional constant. 
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In Section H] we stipulate on the role of A as an accumulant of the 
energy cascade for weak solutions to the Euler equations. We define 
the energy flux through scales £ q by U g = (n q ) where density 7i q is 
the trilinear term that contains all relevant interactions involved in 
the energy transfer. Due to a localization of Tl q observed in [2] all 
appreciable interactions are frequency-local near the £ q scales. We show 
in Theorem 14.31 that the complement of A takes passive part in the 
cascade process in the sense that there is a nested sequence of sets 
G p —> A, in fact G p = Uk >p A k , such that for all p > 

q— »oo 

In other words, the turbulent cascade tends to dump the energy on 
the carrier A. The results of this section are established under the 
natural Onsager regularity condition e = limsup^^dMgl 3 )/^ < oo 
which is known to be suitable for the turbulent interpretation of the 
field (see J2J El 1231 EH HU HI]). Any better regularity of u, namely 
simply requiring e = 0, erases turbulent energy dissipation. This lack 
of dissipation may also occur locally in a proper sense. Therefore it is 
necessary to further confine the cascade to the smaller set A fl S, where 
S is the Onsager singular set of u, i.e. S is the smallest set on the 
complement of which e = 0. This reduction is performed in Section H~3l 
We will demonstrate by several examples that the regularity-related 
intermittency due to 5* and the cascade-related intermittency due to 
A are essentially independent. Therefore one has to consider both to 
have a more complete description of the accumulant. 

The active regions A q can be used to give a more detailed description 
of a measure-theoretic support of the anomalous energy dissipation in- 
troduced by Duchon and Robert in [10]. Recall that for every weak 
solution u G L 3 (Qt) to the Euler equation one can associate a distri- 
bution V{u) such that 

(5) 8 t (lW 2 J +div (u (^H 2 +P^ +V(u)=0. 

In our notation, and under the assumption e < oo , T>{u) is obtained 
as a weak limit of the sequence {TT q } q which is uniformly bounded in 
L 1 (f2y). Thus, T>(u) is a measure of bounded variation. One can easily 
show that suppV(u) C S, but it is unclear whether the topological 
support of T>(u) is contained in the accumulant A or is at all relevant 
to the intermittent cascade. Instead, we will show that the measure- 
theoretic support of T>(u), in the sense of Hahn decomposition, can 
in fact be described by limits of algebraic progressions of the active 
regions A q in a sense defined in Section [61 However, due to a possible 
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slow convergence rate of such a limit, T>(u) may not enjoy the same 
dimensional bound as A itself. 

2. Phenomenology of active volumes 

2.1. Dimensional Fourier analysis on a periodic box. Let B(0, r) 
denote the ball centered at of radius r in IR 3 . We fix a nonnegative 
radial function \ belonging to Cq 3 (B(0, 1)) such that x(0 = 1 f° r 
ICI < 1/2. We define 0(0 = x(£/2) - X (0 and let 

Mx) = f <P(0e 2m( - x dx, 

the usual inverse Fourier transform of on IR 3 . We now introduce a 
dimensional dyadic decomposition. Let L > be fixed length scale 
and let </?_i(£) = xW and = 0(Lf/2«) for g > 0. We have 

E g >-iV^(0 = 1, and 

(6) supp ip q fl supp y? p = for all |p — g| > 2. 

Let T| denote the torus of side length L. Let us denote by J 7 the 
Fourier transform on T|: 

F{u){k) = / w^e" 2 ™^ dx, fc G |z 3 , 

then 

^(/x*) = E /We 2 *** 

determines the inverse transform. We define 

Here ip q must be understood as restricted on the lattice jl? '. In partic- 
ular, U-i(x) = jjs f T3 u(y)dy. It is informative to note the formula for 

h q = T~ x tp q , q > 0, which can be derived via the standard transference 
argument: 

(7) hq(x) = 2^J2 h ( 2 -T + 2<ln )> x ^ J l 
Then for all r G [1, oo] we have 

(8) ||Mr= [pj^ \h q (x)\ r dx^j -2^-^ 
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which is independent of the scale. From (jSJ) follow the classical Bern- 
stein's inequalities, which are also scaling invariant: 

(9) ||u ff || r » < c2 3<? (^~^ 7 )||wj r , 

for all 1 < r' < r" < oo; and the differential Bernsteins's inequalities: 

(10) ||^«,|| r ~(2VL)^||tt,|| r , g>0. 

We will use the special notation X q = 2 q /L for the dimensional wave- 
numbers in order to distinguish them from the a-dimensional ones 2 q . 
The corresponding dyadic length scale is denoted l q = A^ 1 . 

2.2. Active volumes. We assume that a turbulent fluid fills a periodic 
box of linear dimension L. Let us denote it by T|. The motion of the 
fluid is driven by an external stirring force /, which we assume to be 
time independent and of scale rjf ~ L. More specifically, supp J-{f) C 
B(0,c/L) for some a-dimensional c > 1. We observe the fluid on a 
time interval [0, T] long enough to capture the necessary statistics. Let 
us envision that the turbulent motion of the fluid at scale X q consists 
of actively interacting eddies and these eddies fill a region of volume 
V q , which we also call active. Let us now proceed using the classical 
phenomeno logical argument, giving every concept its Littlewood-Paley 
analogue. Our immediate goal will be to derive the following explicit 
formula for V q : 

Here the bracket (•) denote the average in space-time on the domain 
Qt = [0, T] x T|. Let us start by noticing that since tp q has mean 
zero, for q > 0, we can prescribe the turn-over velocity of an £ q -eddy 
at location x to be u q (x). Let U q be a characteristic velocity of an 
£ q -eddj. In order to single out active eddies from passive ones, one 
will use an LP- average such as U q ~ (\u q \ p ) 1 ^ p . The minimal value of 
p proved to be suitable for studying intermittent cascade turns out to 
be 3 (see Section H]), while higher values can be adopted to formulate 
multi-fractal hypotheses as shown in Section 14.51 So, we define 

(12) u q = -^(Ki 3 } 1/3 , 



Vq 



or more explicitly, 



(13) U q = [ — I \u q (x,t)\ s dxdt 



1/3 

3. 
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Even though the integration in ( |T3l) is performed over the entire do- 
main, we will prove later that there exists indeed an active region 
A q C VL T with \A q \< TV q} and 

(14) / \u q \ 3 dxdt ~ / \u q \ 3 dxdt. 

J Aq JQt 

The input energy produced by the force / is passing from larger to 
smaller scales. The energy flux per unit volume carried through the 
scales of order £ q is given by 

(15) ^ = ^ = ^<W S >- 

**q ^q v q 

On the other hand, 

K. 



(16) e 
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where K q is the average kinetic energy of an active £ g -eddy given by 
(17) ^ = fW), 

Vq 

(here again we use the active proportion of the volume) and t q is the 
typical turnover time given by 



<l ^qVq 



1/3 



(18) tq ~ u q ~ L{\u g \*y/z- 

Putting together (TTTj) and ( 1T81) we obtain another expression for e q : 

(19) ^ = -|^(kl 2 )(KI 3 ) 1/3 - 

"q Vq 

Equating f lT5|) and f|19|) we finally obtain ( TTTj) . 

2.3. Intermittent dimension. We define the co-dimension f) as an 
asymptotic scaling exponent: 



(20) () = liminf 



log 2 (L 3 /^) 



q^-oo q 

Let the dimension be d = 3 — 1). Observe that since all volumes are 
bounded by L 3 , we have d < 3. The dimension <i is defined without 
any scaling hypothesis like in the /3-model and therefore is a property 
of the field. 
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3. Active regions 

Let us assume as before that u(x,t) is the velocity field of a turbulent 
fluid in the periodic box T| on a time interval [0, T], and V^'s are defined 
by ( II ip . In this section we will introduce regions that capture most of 
the active eddies, i.e. eddies with at least a certain turn-over speed s q . 
In order to properly define the speed and the regions we will make use 
of atomic decompositions. Let us briefly recall the definition (for more, 
see PP). For q > and k G [0, 2 q — l] 3 D Z 3 we define the dyadic cubes 

Q qk = [0,Q 3 + £ q kcT 3 L , 

and dilated cubes 

Q*qk = (J Qqk'- 
|fc'-fc|<3 

For any M > 1 one can find a spacial decomposition of u(t, •) into 
M-atoms 

oo 

(21) u(t, x) = } j s qk {t)a qk {t,x), 

g=o fce[o,29-i] 3 nz 3 

where supp(a g fc(t, ■)) C Q* fe compactly, and 

max {A- |/9| ||afa 9fc ||ioc (nT) } < 1. 

For every r < oo one has 

(22) £ M*)rl ~IM*,0llr. 

Let us recall how the coefficients s qk are constructed. Let Xqk be the 
characteristic function of Q qk and let rj be an approximative kernel 
supported on 5(0, 1). Let r] q (x) = 2 3q rj(\ q x), x' g k = X q k*Vq and b qk = 
u qX' qk - Then let 

s qk (t) = max{A-^||9fV(t,-)||oo} 

<^gfc bqk I Sqk- 

Let us now fix any positive a-dimensional decreasing sequence a q — > 
with zero exponential type 

lim ^ = 0. 

g^oo q 
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Definition 3.1. A active region occupied by eddies of size i v at time 
t is defined by 



(23) 

Furthermore, let 
(24) 



A q (t) 



u 



Qqk- 



k:\s q k{i)\>cr q 



<l"ql 3 > 



A q = {(t,x) : x e A q (t)}. 



\A q \ < ca-'%T, 



\u q \ dxdt < 



\u\ dxdt < 



\u q \ dxdt, 



In the following theorem we collect several properties of active re- 
gions that are independent of particular nature of the flow or even the 
evolution law of the field u. For a set A C Qt, let A(t) = {x G T| : 
(x,t) e A}. 

Theorem 3.2. Let A q be the active region defined by ([231 . and let 

(25) A = lim sup A q = f^ =1 U q>p A q . 

q—>oo 

Then for some absolute constant c > we have 

(26) 1 A 1 ^ —" 3l 

(27) 

(28) dim n A(t) < d, for a.e. te [0,T]. 

In particular, if d < then A(t) = for a.e. t. 

We thus see that the active regions are regions where most of the 
characteristic velocity U q is concentrated, yet the measure of A q on 
average does not exceed V q up to an algebraic multiple. Inequality ( 1281) 
confirms the dimensional prediction of the /3-model, and we will argue 
in the next section that A indeed represents the set where the energy 
flux accumulates. 

Proof. Let N q (t) be the number of cubes in the union (|23|) . Then 
trivially 

27 



\A q (t)\<-N q (t). 



On the other hand, 



*3 



Nq(t) 



u 



Q 



qk 



k:\s qk (t)\>^^ 
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Integrating in time we obtain 

(29) 

In view of ( )22l) we have 

1 r T ^ cr 3 (In l 3 \ 3 r T 

/ M'dxdt- £M*)| 3 ^>^)^/ N q (t)dt. 

>n T \ J0 l \ \\ u q\ I JO 



- jf N q (t)dt < \A q \ <-y o N q (t)dt. 



Thus, in view of 



as claimed. 

Now let us observe that on the complement A c q = Qt\A, we have 

\u\ 



/I 1 3 \ 

< ca ) , y . We therefore obtain 

— 9 (KI 2 ) 



/" /" in |w| 3 cfe<it r 

/ ix cfrrdi = / it \u\dxdt < a q r T — — — — / ix dxdt 
J a- J A c J Qt \u\ 2 dxdt J A c 



< o„ I |ix| 3 dxdt, 

'tor 



which proves ([2] 

Let us first recall that the Hausdorff dimension of a set A C T 3 , 
dim^(A) is the smallest d for which the Hausdorff measure vanishes 

H d+S (A) =limH d+s , e (A) = 0, 



£->0 

for all 5 > 0, where 



"H dj£ (A) = inf <^ ^(diam Ai) d : A C U^, diam A, < £ > . 



i=l 



Directly from ff26|) we have 



(30) ^j o T N q (t)dt<a; 3 X' q V, 



Given 5 > 0, we note that 



(31) \ < 2-^- 8 ' 2 \ 

for all q large enough. Let us now assume that d > 0. Then 
H d+SA (A(t)) < n d+5/q (U p>q A p (t)) <Y,N p (t)j^. 

p>q P 
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Integrating in time, using ( )30l) . ( 13TT) and the fact that a q — > alge- 
braically, we obtain 



f n d+SA (A(t))dt < ^5 < T E 

JO ~Z, °nK ~Z, 



L d 2 8/2 



p>q "p-^P^ p>q a P^P 

as q — > oo. So, in the limit we obtain 

H d+s {A{t))dt = 0. 

Hence, dim^ A(t) < d + 5 for a.e. i e [0, T], which concludes the proof. 
If d < 0, then and (EU) imply 

-3o-g<5/2 



/ ^iV p (t)rft<a- 3 2 
- 70 P>q 



for q large enough. Thus the measure 



P>Q J 



and on the complement the set A q (t) is empty. Passing to the limit we 
conclude that A(t) is empty a.e. □ 



4. Energy dissipation set 
Let us consider the Euler equations in the periodic box T 3 L : 

(32) — + (u • V)« = -Vp + /, 

(33) V • u = 0. 

A vector field u G C w ([0, T}; L 2 (Tf )), (the space of weakly continuous 
functions), is a weak solution of the Euler equations with initial data 
u e L 2 (Jl) if for every i/; e Cg°([0,T] x T|) with V x • ^ = and 
< t < T, we have 
(34) 



u-ip— UQ-tp — u • d s ip = / (u ® u) : Vip + f ■ if), 

t|x{*} Ji 3 L x{o} Jn t J n t Jn t 

and V x • u(t) = in the sense of distributions. We define the operation 
: by 

A: B = Tv[AB}. 
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4.1. Energy flux and density. In order to properly define the flux of 
kinetic energy across the scales, let us fix a q > 1 and test against 
the filtered field 



9-1 



«< 9 



P =-i 



Denoting w <(? = J r_1 (^/x(2 9 • /L)F(u)) which also represents a dyadic 
filtration, we obtain the following budget relation for the energy dissi- 
pation rate across the wavenumber \ q : 



(35) 



I (u (g) u) : Vu <q dx + f ■ udx, 
Jt 3 t Jn 



(notice that the latter integral is independent of q for q larger than the 
integral scale). Let us denote 



T\- q (t) = (u®u): Vu <q dx 



the total energy flux due to nonlinearity. It represents the averaged 
contribution of all nonlocal interactions to the energy exchange across 
the wavenumber X p . By antisymmetry of the trilinear term, we have 

/ {u® u <q ) : Vu <q dx = 0, 
M 

which is a statement of the fact that eddies larger than £ q on average 
do not carry the energy across the scale l q . Thus, 

(36) U q (t) = / (u®u> q ):Vu <q dx. 

Jn 



L 



In this form the flux is clear from the large scale " shuffling" . We notice 
that the Fourier support considerations reduce the formula f)36p further 
to 

(37) n ? (£) = / (u p i <S> Up") : Vupwdx. 

L p'>?-i, p">g, p"'<q 

\p'-p"\<2 

This formula shows detailed contribution of scales to the energy budget 
relation. We will see in Lemma 14.11 that the remote scales do not in 
fact have a considerable impact on the flux due localization (see [2] and 
Lemma [4. II below) . We proceed now to define the energy flux density 
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as the integrand of ( 1371) : 

(38) n q = {u,pi ®u p n) : Vu p m. 

v'>q-l, p">q, p"'<g 

\p'-p"\<2 

Thus, 

(39) U q (t) = / TT g (x,t)dx. 

hi 

In addition we define for each K < q the truncated density 

(40) tt^ = ^ (u p i <S>u p ») : Viy». 

p',p",p'" as above 
p',p",p'"e[q-K,q+K} 

Following the ideas set forward in [2] we obtain a localization prop- 
erty of the flux in the next lemma. 

Lemma 4.1. Suppose that 

(41) e = limsup X q (\u q \ 3 ) < oo. 

q— >oo 

Then we have 



(42) limsup(|7r 9 |) < ce, 

/or some a- dimensional c > 0. Moreover, 



(43) lim limsup(|7r, - 7rf |) = 0. 

Proof. Using the Holder and differential Bernstein's inequalities, we 
have 

<W> < E (lvl 3 ) 1/3 (lvl 3 ) 1/3 V(lv"l 3 ) 1/3 



p'><?-L p"><?, p"'<9 

|p'-p"|<2 



<(e<m 3 >) (e*j 

\P>9-2 / \p<g 



2/3 / x 1/3 

3\ 



l«p| d > 



2/3 / x 1/3 

£ Wvn 3 >) E A ?V A P (i^i 3 ) 



\p>q—2 / \P<q 



Clearly, the sums above are of convolution type with the exponen- 
tially decaying kernels A^A" 1 and A 2 A~ 2 . Hence, f|42|) follows. If one 
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of p',p",p"' is outside [q — K,q + K], then in at least one of the sums 
above the summation is performed over p's with \p — q\ > K. Thus, 

limsup(|7r g - 7r* |) < ceA A , 2/3 , 

q— ¥00 

and (j33D follows. □ 

4.2. Energy flux and active regions. In this subsection we will 
connect the energy dissipation with active regions. 

Definition 4.2. We say that the energy flux concentrates on a mea- 
surable set G C Qt, or G is an energy dissipation set, if there exists a 
decreasing sequence of open sets G q such that G = fl^G^ and 

( 44 ) }™^ P \n T \ G J) = 
for all q G N. 

We will now prove that the energy flux concentrates on the active 
region A defined in Theorem 13.21 

Theorem 4.3. Suppose u G C w ([0, T}; L 2 (T|)) is a weak solution to 
the Euler equations, and 

(45) e = limsup X q (\u q \ 3 ) < oo. 

q— >oo 

Then the energy flux concentrates on A = limsup A q , where A q are 
the active regions. 



Proof. Let G q = U^ =q A p and let G c q = £l T \G q . According to LenrnaE 
it suffices to show the limit 

(46) hmdvfl |) = 0, 

p— s-oo i 

for all K, q G N. For all p > q + K we obtain the following estimate 

(Wp \ G o\) < Yl (\ u p'\go\ K"\ G o\ iVvlogl) 

P ',p",p'"e[p-K,p+K] 

< ^ ( I V U«, I I U P" I A=„ I I V V" I A= /w I ) • 

' * p' p" p'" 

p',p",p'"e\p-K,p+K] 
Let p* G {p',p",p"'} be such that 

(|wp»| 2 ) = max{(|M p /| 2 ), (|vl 2 )> (IV'I 2 )}- 



/// 



We will consider three different cases: p* = p', p* = p", and p* = p 
The first two are similar, so we consider only the first and the last ones. 
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Assume p* = p' . We have 



(Ivl 3 ) 



\Up>(t,x)\ A c | < co v ,- 

p' {\Up> 



Then we obtain 

(lvU c , 



Up" I fic. 



Nil 



P"'\A C 



,\)<*P 



(M 3 ) 
(M 2 > 



(\up»\ |Vty»|) 



(l M p'| 3 ) , it |2\ 1/2/1 , 1 2\ 1/2 

~ a p'7i — mV'dvl ) (Ivl ) 
(Ivr) 

< 0~pi\2KXp'(\Up>\ 3 ) 

as p — > oo. 

In the case where p* = p'" we notice that for every (t, x) G A c m , 



|Viy//(t, x)| < 



< \pniUpi 



22 Sp»'k{tNa p »i k (t,x 
k 

So the argument above works as well. We thus obtain 

P 



(IV"| 3 ) 

(Ivl 2 ) 



K\ G o\) <c{K)ea p _ K ^0, 



as p — > oo. 



□ 



We note that Theorem 14.31 provides no guarantee that the energy 
dissipation will occur if A ^ 0. In fact at present no example of 
a weak solution to the Euler equation is known to possess a non- 
trivial flux (7f q ) -/+ 0, and at the same time to have Onsager regularity 
sup„ \ q (\u q \ 3 ) < oo. In the next section and Section [5] we will find 
active regions explicitly for some known singular solutions. 

4.3. A further confinement of the energy flux. Let us consider 
the following two dimensional example. Let e} be the vectors of the 
standard unit basis. We fix a large s > and define 



u(x 1 ,x 2 ) = ^2 



u 



<?=! 



U, 



,, - e 2 — sin(A 5 Xi). 



Clearly u is a smooth stationary parallel shear flow. By a direct com- 
putation for this example we have V q ~ 1, so d = 2, s q j. ~ A^ s and 
hence, A q ~ Qt, A = Qt ■ Yet, it is clear that the energy dissipation 
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set is empty. The reason for this is that the conclusion of Theorem 14.31 
is based solely on the intermittency character of the flow, i.e. satu- 
ration of the L 3 -average relative to the L 2 -average. This particular 
example lacks it completely - it is uniformly "active" throughout the 
flow domain. So, in order to make the prediction of the theorem more 
precise we need to dispose of the regions where regularity of u exceeds 
Onsager's threshold, namely e(u) = limsup^^ X q (\u q \ 3 ) = 0. In order 
to state the local version of e(u) we use the terminology introduced in 

E21. 



Definition 4.4. Let u : [0,T] — > L 2 (T 3 L ) be a time dependent field. 
We say that u is (Onsager-)regular if e{u) = 0. Now let U C Qt be 
a relatively open set. We say that u is regular on U if (mp) is regular 
for every scalar function ip G C£°(U). We denote by D the union of all 
open U's on which u is regular, and call S = Qt\D the singular set of 
u. Clearly, S is closed. 

As shown in [22] for every weak solution to the Euler equation u G 
C w ([0, T]; L 2 (Y\)) the local energy equality holds on its regular set D: 



j \u\ 2 (p- / MV- / M 2 <9^ = / (\u\ 2 + 2p)u ■ Vy?, 

JD(t") J D(t') JD JD 



(47) 

lD(t") JD(t'] 

for all (p G C^°(D) and t',t" G [0, T\. In terms of anomalous dissipa- 
tion distribution introduced by Duchon and Robert [10] it means that 
suppP('u) C S. 

Theorem 4.5. Under the assumptions of Theorem \^.3\ the flux con- 
centrates on An S . 

Proof. Let S q = {y G Vl T : dist{y, S} < ±}. Then S = f) q S qi A n S = 
r\ q (S q n G q ), and so, it suffices to show that 

(48) lim(K| |) = 0, 

for all K, q G N. Following the proof of Theorem 14.31 we argue 

(Irf |g§us«I) ^ (IvIgcI \ u p"\gc\ |VV"|g§I) 

p',p",p'"e\p-K, P +K] 

+ Yl (Msfl IvUjI l v v|5=i)- 

p'yy"e[p-x, P +x] 

It thus suffices to prove convergence of the second sum to zero. Since S° 
is closed and disjoint from S we can find a cut-off function (p G Cq°(D) 
such that = 1 on 5° We use the following representation 

(49) (uip) p = Upip + uip p + r p (u,tp), 
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where 

r p (u,<f)(x,t) = / h p (y)(u(x-y,t)-u(x,t))(<f(x-y,t)-<f(x,t))dy. 

M 

Restricting ()49l) onto S° we obtain 

u p '\ S c = (uip) p >\ Sc - uip p >\ S c - rp/(u,ip)\ S e- 
We then trivially estimate 

(|u p /| 5 c| |w P "|5c| |Vtt p '»| 5 c|) < (\{mp) p t\ \upii\ \Vup"i\) 

+ (\u(p p '\ | u p ii | |Vw p "'|) 
+ (|ry(u, <f)\ \u p »\ \Vu p m\). 

For the first term we have 

(|(^VI |vl IVV'I) < C(K)(A p ,(|(^V| 3 )) 1/3 (A p »(|vl 3 )) 1/3 x 

x(A p /»(|v'| 3 )) 1/3 > 

which vanishes as p — > oo by the regularity of zap. The second term 
vanishes trivially by the regularity of ip. As to the third, we have by 
Minkowskii's inequality 

(IVK^)I 3 ) < / \h p (y)\(\u(--y,t)-u(-,t)\ 3 \<p(--y,t)-<p(-,t)\ 3 )dy 

<c(|u| 3 ) || V^loo / \h v {y)\\y\dy<C(ui<p)\-\ 
Using this estimate we obtain 

(Mu,<p)\ ivl |Vvl> < c7( m ,^)a;, 1/3 ((|vI 3 )) 1/3 (a P '»(Iv"I 3 )) 1/3 

<C(u,v,K,e)\;, 1/3 ->0. 
This concludes the proof of the theorem. 

□ 

4.4. Application to the scaling laws of turbulence. In this sec- 
tion we derive some bounds that replicate the classical power laws of 
turbulence with intermittency correction in terms of 

q 

Note that 2 qdq represents the number of active eddies of size l q . More 
precisely, thanks to ( 1301) we have 

i ^ N q {t)dt<o~\L\) d «. 
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Like before, < stands for an inequality up to a universal constant. 
Let us denote e = sup X q (\u q \ 3 ), and e = inf(|7r 9 |). Then 



i<(W> 



< 



Pll M pll3' 



for all q. Here the second inequality follows from the local estimates 
on the flux (see j2] or Subsection 14 . 3 f) : 

(50) (k|> <J2 K *-p X 

v 

where K is the following localization kernel 

(51) K q = 



\ 2 q '\ q < 0; 



A," 473 , g>0. 

Now let us recall that the energy spectrum in the Littlewood-Paley 
settings is defined as follows (see [6]) 

(52) E, = M>. 

To take into account intermittency we notice that 

(53) V q = L d "X d q "- 3 
and obtain 

<M 2 > = (\uf) 2/ X /3 L- 1 = <K| 3 > 2/3 (^)V3-i < L E 

Xq 

Hence we recover the energy power law with intermittency correction: 

(54) E q < 



^2/3 



xT(Lx q y- d ^' 

Note that this bound recovers ([1]) with d = sup d q . 

Let us now estimate the energy spectrum from bellow. For this we 
use the estimates on the flux fl50|) obtaining 

£ 2/3 <<m 2 / 3 < (e^-a(ki 3 )) 1 <E*f(i^f) 2/3 

\ V J p 

= E<-lK /3 (\u P \ 2 K /3L = Y,<-lK /3 ( L K) 1 - dp/3E ^ 

V V 

Combining with (I54p we have the following two-sided estimate: 
(55) < * {Af {LX q f- d ^E q } q < e 2 '\ 
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Notice that since the kernel K 2 ! 3 has an absolute weight concentrated 
mostly near 0, the lower bound can be interpreted as the averaged spec- 
trum over nearby wavenumbers. Such an averaging is not uncommon 
in experimental science. 

Next, let us consider the velocity displacement S y u(x,t) = u(x + 
y, t) — u(x, t) and for < £ < L define the generalized isotropic second 
order structure function by 

S 2 (£) = -H- f (\Sieu\ 2 )d6. 

Lemma 4.6. Let e < oo and d — 3 — f) be defined as before. For every 
5 > there exists an adimensional C$ > such that 



2 „2 



(56) S 2 (£) < es£s 



1 d z 4 

^a" 5 /P\3 



l, +c Hz 



for all £ < L. 

Proof. Let us fix y ^ and using (|53|) observe for all q > 

(IM 2 )<EM 2A ^KI 2 ) + E(Ki 2 ) 

p<q p>q 

(57) s ^ E \vKK /3 K 2/3 + ^ E v i" x 



2/3 

/' 

p<q p>q 



< ^ E IvW^)*- 1 + e 2 / 3 E( LA p)^ _1A p 2/3 

p<q p>q 

Let 5 > be fixed. Then in view of lim sup^^ d q = d, there exists a 
p G N such that for all p > p one has d p < d + 5. Let us suppose that 
\y\ < j—. Then we choose q = [\og 2 (L\y\~ 1 )} + 1 > p . Notice that 

X q ~ l/|y|- Then splitting the first sum in (j57p we continue 

(58) 

(\^u\ 2 )<C 5 e 2 / 3 \y\ 2 L^ 3 + e 2 l 3 \y\ 2 E (LA,)"/ 3 "^ 

P0<P<<? 

+ e 2 / 3 E(^A p ) d/3 - 1+5 A p 2 / 3 

< C 5 e 2 / 3 |2/| 2 Z.- 4 / 3 + £ 2/3 M 2 Af (L\ q f s - 1+S + e 2 / 3 (L\ q ) d / 3 - 1+s \ 

< C^yfL-W + s 2 /\LX q ) d ^ +s X7 2 /\\y\ 2 X 2 q + 1), 
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where C$ > is a constant dependend only on 5 > 0. Now recalling 
that X q ~ l/\y\, we obtain 



<M 2 }<? 2/3 M 2/3 



4/3 



l-d/3-<5' 



which is dSB]). 

Now let us suppose that \y\ > Then we choose q = p and 



continue from ()5j 

(59) (\S y u\ 2 ) < C & e^\y\ 2 L^ + e^(L\ Po )l-^\-^ 

Notice that the power of A po is negative, so using our assumption we 
replace it with 



(\S y u\ 2 )<C s e^\y\'L-^ + e^ O^j 



l-d/3-5 



£ 



2 /3 1 y 1 2/3 



4/3 



1— d/3— 5 



2/3 



which immediately implies ([56 



□ 



By a similar argument using the fact that e < oo we find (l^wl 3 ) < 
Thus, by interpolation, for any 2 < p < 3 we have up to a higher 
order correction term 



(60) 



s P (e)<ek^, c P = | + (3 - rf) (i - 1) • 



Formula flBTJ]) is exactly the one that appears from the dimensional 
argument of the /3-model. We note however that it is not consistent 
with experimental results for larger values of p. A finer tuning can be 
achieved by using the multi- fractal formalism of Frisch and Parisi [21]. 



4.5. Remarks on multi-fractality: future research. As we have 
seen, no uni-scaling hypothesis used in the original /3-model [2] is 
necessary to justify some intermittency corrections for deterministic 
fields. In fact, the accumulant A can be used as a base to build a 
multi-fractal theory similar to the one proposed in the celebrated work 
of Frisch and Parisi [21] . In our context this can be achieved by defining 
a nested sequence of active volumes, regions, dimensions, etc, based on 
saturation of the L 2 -L p averages. Indeed, if the argument of Section I2~2l 
is started by using //-averages of the velocity field as characteristic 
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speeds, one obtains 

(61) VV=L* M * . 

(\u q \*)p-> 

The threshold speed for an active eddie becomes 

(62) |4 



.../■>, „<K!!>f: 

(\u q \ 2 )p- 



One can define active regions Aq and accumulant A^ in a similar 
way. By interpolation, one has A q p " ) C A q p ' ] for all p" > p' > 3, and 
hence {A^}s<p<ca defines a foliation of our base set A. Multi-fractality 
can therefore be viewed as a set of scaling parameters: the dimensions 
as co-exponents of the volume sequences Vq P \ and a spectrum of 
scaling exponents of the field on accumulants A^ p \ 

In a subsequent work the analysis above will be carried out formally, 
but the basic conclusion remains the same: multi-fractality, just like 
the other basic concepts introduced in this note, is not an assumption 
imposed on a field - it is a property of the field determined by its 
volumetric quantities. 

5. Examples 

In this section we examine several examples of fields to illustrate 
how the active regions can be computed explicitly. We also show that 
generally there in no inclusion of the carriers A and S in either side. 

Let us consider first the stationary vortex sheet solution, where 
u(x,y,z) = (H(z), 0, 0), where H is the Heaviside function. Then 
by a direct computation, (\u q \ p ) ~ £ g , for any p > 1 and q > 1. Thus, 
e G (0, oo), V q ~ £ q , and therefore d = 2, which is exactly the dimen- 
sion of the set of discontinuities for u. Notice also that S = {z = 0}. 
Furthermore, A = {z = 0} as well. One can see it by manufactur- 
ing atoms for u out of smoothed characteristic functions of the dyadic 
cubes. Then s q k become comparable to 1 only about the 1/A 9 vicinity 
of S, and zero away from it. This makes A q be a sequence of l/A^-thin 
slabs converging to S. Thus, in the limsup A q — > S. Note that u is sta- 
tionary, / = 0, and therefore the energy conservation holds, despite of 
the fact that neither A or S sees it. Even more generally, one can show 
that for any smooth vortex sheet A and S coincide with the sheet, yet 
the energy conservation holds. This is due to the particular kinematic 
condition on the sheet (the particles cannot cross the surface of the 
sheet) that follows from the weak formulation of the Euler equations 
(see [22J). 
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Computation of the active regions in the case of a one-point sin- 
gularity can be done explicitly too. We recall for a moment a two- 
dimensional example studied in [23] where u is assumed to be —1/3- 
homogeneous in the radial direction near the origin, it is a stationary 
solution to ( 132|) with smooth forcing /, and thus u = V~ L ?/> where 
ip(r } 6) = r 2 / 3 fy(9) near the origin, far from it, and ^ G C 2 ([0, 27r]). 
Here (r, 9) stand for the polar coordinates. Clearly S C {0} and it is 
easy to verify directly by integration that (\u q \ p ) ~ i q 2 p ^ 3 \ for p < 3. 
Thus, S = {0} and the solution overall is Onsager-critical. Further- 
more, V q ~ Eg, and d = 2 — t) = as it is supposed to be. To find 

i i 1/3 

A we note that the condition \s q k\ > cr q X q can only be satisfied near 
the origin since away from it, s q k — > exponentially fast in q. Thus, 
A q — > {0} = A. The smoothness of / necessitates the singular part of 
the pressure to be given by p = r~ 2 / 3 P, P is a constant and ^ satisfies 
the following ODE 



As shown in (23], (1631) implies II = 0, and hence all such stationary 
solutions have no anomalous dissipation, again in spite of the fact that 
A = S 7^ 0. It is not clear at the moment whether there exist at 
all Onsager-critical weak solutions to the Euler equation with smooth 
forcing, and a one-point singularity that have non-trivial anomalous 
flux. 

The analogous three dimensional construction can be considered as 
well, where in spherical coordinates u = -^j3U(9,(j)) near the origin, 

and U is smooth vector field. In this case we have (\u q \ p ) ~ if 2p ^\ 
for p < 3, and hence V q ~ £ 3 confirming again the dimension d = 
is consistent with the structure of the singularity. At this moment an 
examination of the flux as in 2D has not been performed although we 
have reasons to believe that the result remains the same. We note 
that the analogue of the ODE ( 163]) becomes now a system of coupled 
nonlinear second order PDEs on components of U. 

So far we have seen examples with the inclusion S C A. Generally, 
this is not the case. Let us consider cubes of integers C q of side width 



(63) 



3(^') 2 + 4 ^ 2 + 6W = P. 



The flux is given by 



(64) 
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Ag/10 placed at the frequency \ q e\, and define 

i 

T{u q ){k) = %-— XCq (k). 

Aq 

Then ||w 9 ||3 ~ l/A^ 3 , ||w g ||2 ~ 1/Ag^ 6 , and hence V q ~ 1/A 3 , and d = 0. 

2/3 

The threshold for active eddies becomes \s q k\ > <J q \ q , yet one can 
see that u q is the product of 3 modulated and scaled Dirichlet kernels 
\ q 7 ^ 3 Df 3 ^ 10 . Thus \u q \ ~ Xq at the origin and u q concentrates near 
the origin as q — > 0. So, by choosing a sequence a q decaying slowly 
enough we can conclude that the eddies over the threshold concentrate 
near the origin as well. This shows that A = {0}. Now let U C T| 
be open and if G C^°(U) and if > 0. Let a = J > 0. We have 
(u(p) q ~ MgV 2 by the same analysis as in the proof of Theorem 14.51 
Then due to the smoothness of ip and uniformity of Fourier modes of 
u q in the box C q we obtain 

(65) F(u q <p) (k) = F{u q ) * F{<p) ~ a/ Aj/ 3 

for g large enough and for all A; in a sub-box of C q with edges of size 
A 5 /10 — N. The larger the N the more precise the identity (|65|) becomes 
in the limit g — >• oo. Thus, up to a negligible error, 

ll^lls > HA-^fVivlls ^ «Af ■ 

This shows that uip is Onsager-singular, and thus S — T|. 

Finally we note that an example is available with smoothness 1/3 in 
the L 18 / n -average sense (u G B\^j n ^ in the notation of Besov spaces) 
that solves the Euler equation with smooth forcing / = (0,0,cos(x)) 
and has anomalous dissipation (/ • u) ^ 0. The details of this construc- 
tion will be presented elsewhere. 

6. Connection with Duchon-Robert's approach 

Let u be a weak solution to the Euler equation with e < oo, and 
let T>(u) be the distribution of Duchon and Robert as described in the 
Introduction. From Lemma 14.11 we can see that {7r g } g is uniformly 
bounded in L 1 (f2y). By following the same calculations as in [10] one 
can show that in our terms T>(u) is the weak*-limit of this sequence, 
hence T>(u) is a measure of bounded variation. In this section we will 
use the active regions to describe a measure-theoretic support of V{u) 
in the sense of Hahn decomposition. 
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We first make several definitions. Let B C £It be a measurable set, 
and fc > be an integer. We say that B is a fc-cluster of the sequence 
if for every open set U containing B one has 

(66) lim max{|gi - q 2 \ : q. L > Q, (U qi < q < q2 A q ) D U = 0} = k. 

In other words for very U the gap in the string of indices q such that 
A q n U = eventually does not exceed k, and fc is the minimal such 
gap. If k — oo, we define a oo-cluster as a set B which is not a fc-cluster 
for any finite k > 0. 

Let B be a closed fc-cluster. Denote k = k(B). Clearly, if B\ C -B2, 
then k(B\) > k{B-i). We call 5 minimal if for every proper closed 
subset .B' of -B, one has k(B') > k(B). By Zorn's lemma every closed fc- 
cluster set B, with finite k, contains a minimal fc-cluster subset. Indeed, 
let {B a } a( zi be a chain of fc-clusters ordered by inclusion. Let B = 
^aeiB a - Then for every U containing B , B a C U eventually, i.e. for 
all a > a . Since B a is a fc-cluster, (1661) holds. This shows that S is a 
lower bound for the chain, and Zorn's lemma applies. 

Consider the set 

M= |J |J B. 

0<fe<oo B is minimal fc-cluster 

Theorem 6.1. Suppose u is a weak solutions to the Euler equation 
with e < 00. Then 

(i) Q T \M is a null-set ofV(u); 

(ii) supp£>(w) C M C n g >iU p>9 A p . 

Proof. For (i) it is enough to show that T>(F) = for every closed 
subset of Qt\M. Since FC\B = for any minimal fc-cluster set B, and 
F is closed, F itself is not a fc-cluster for any finite fc. Hence, for every 
fc > there exists an open U D F and for any N > there exists a 
string g^, > iV with \q' N — q'^\ > 2fc + 1 such that A q (~)U = 0, 
for all q' N < q < q%. Let gjy = [(^ + Qn)/ 2 } and let G C^{U) be 
arbitrary. Then 

/ (p7r qN dxdt — > / 

as —7- 00. On the other hand, by the same computations as in the 
proof of Theorem 14.31 and since there are no active regions in U with 
indices near q^, we obtain 

< Co q , 0, 



7T 



'In 



(fidxdt 
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as 



iV — > oo. Thus, in view of Lemma |4JJ 



<pdV 



< oil 



where o(l) -»■ as fc -> oo. This implies |X>|(F) < |X>|(C/) < o(l). 
Since F is independent of fc, letting k — > oo we obtain the desired 
result. 

The conclusions of (ii) are even more straightforward as seen from 
the argument above - there is no clustering of active regions at all 
outside of M. 

□ 
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